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ABSTRACT 

We present a model of a freely precessing neutron star which is then compared 
against pulsar observations. The aim is to draw conclusions regarding the structure of 
the star, and test theoretical ideas of crust-core coupling and superfluidity. We argue 
that, on theoretical grounds, it is likely that the core neutron superfluid does not 
participate in the free precession of the crust. We apply our model to the handful of 
proposed observations of free precession that have appeared in the literature. Assuming 
crust-only precession, we find that all but one of the observations are consistent with 
there not being any pinned crustal superfluid at all; the maximum amount of pinned 
superfluid consistent with the observations is about 10~^° of the total stellar moment 
of inertia. However, the observations do not rule out the possibility that the crust and 
neutron superfluid core precess as a single unit. In this case the maximum amount of 
pinned superfluid consistent with the observations is about 10~* of the total stellar 
moment of inertia. Both of these values are many orders of magnitude less than the 
10~^ value predicted by many theories of pulsar glitches. We conclude that superfluid 
pinning, at least as it affects free precession, needs to be reconsidered. 
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1 INTRODUCTION 

This purpose of this study is to collect the important factors 
that influence neutron star free precession, combine them 
into a single unified model, and then test this model against 
observations. In particular, we wish to learn about the geo- 
logical history of the star, decide whether or not superfluid 
pinning in the inner crust is required to flt the data, and ask 
whether it is the just the crust that participates in the free 
precession, or whether some (or all) of the interior neutron 
superfluid participates too. By 'observations' we refer to the 
handful of pulsars where a smooth modulation in the pulse 
timing and/or structure has been detected. We hope that 
the model and equations presented here will be of use to ob- 
servers when new free precession candidates are discovered. 

Clearly, we can only hope to extract useful informa- 
tion from the observations if our model of neutron star free 
precession is sufficiently realistic. To this end, our model 
includes the effects of crustal elasticity, inertial coupling 
(where the interior fiuid pushes on the surrounding crust), 
and a superfluid component pinned to the inner crust. We 
also examine dissipative processes that tend to enforce coro- 
tation between the crust and core. We will argue that cur- 
rent ideas concerning crust-core coupling suggest that only 
the crust (and perhaps the charged plasma in the fluid core) 



undergoes free precession, with the interior neutron fluid 
simply rotating about an axis fixed in space. This has bee n 
assumed implicitly in earlier works ( Alpar fc Sauls 1988 ), 
but not commented on explicitly. 

The structure of this paper is as follows. In section ^ 
we describe how strains in a neutron star crust can pro- 
vide a non-spherical contribution to the moment of inertia 
tensor. In section |^ we describe our free precession model. 
Dissipative processes which convert the precessional energy 
into heat or radiation are considered in section ^. In section 
^ we describe how the finite crustal breaking strain places 
limits on the shape and wobble angle of a freely precessing 
star. The way in which the signal from a pulsar is modu- 
lated by free precession is described in section ^. In section 
1^ we combine the results of the previous sections to extract 
as much information as possible from the proposed observa- 
tions of free precession that have appeared in the literature. 
Our conclusions are presented in section |. 



2 DEFORMATIONS OF NEUTRON STARS 

To undergo free precession, a star must be deformed in some 
way, so that its shape and moment of inertia tensor differ 
from that of an unstressed fluid body. In this section we 
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describe how strain in the solid crust can provide such a de- 
formation, using the model of Pines & Shaham (1972), who 
in turn made use of the terrestrial analysis as presented in 
Munk & McDonald (1960). For reasons of tractability a ro- 
tating but non-precessing star is considered. We will not con- 
sider deformations due to magnetic stresses, although these 
could be important for the free precession of very slowly ro- 
tating/very highly magnetised stars (see e.g. Melatos 1999) . 
Note, however, that a pulsar whose deformation was due en- 
tirely to magnetic stresses wo uld no t display any modulation 
in its pulsations — see section 6.1.1. 



Begin by writing down the total energy of the rotating 
star. Let /star denote the moment of inertia of the spheri- 
cal star, i.e. the moment of inertia the star would have if 
it was unstrained and not rotating. When rotating, the mo- 
ment of inertia about the rotation axis must be greater than 
this, and so can be written as /star(l + e). We will refer to 
e as the total oblateness. We imagine the crust to have so- 
lidified from a hot, liquid, state in the geologically distant 
past, leaving it with a 'reference' or zero-strain oblateness 
eo- This parameter will then change only via crust-cracking 
or a gradual plastic creep. The star's energy is a function of 
e and eo according to: 
r2 



E — Esta 



2/star(l+e; 



+ Ae^ + B{ 



(1) 



Here -Bstar denotes the energy the star would have if it was 
spherical. The second term is the kinetic energy. The third 
term is the increase in gravitational potential energy due 
to the star's shape no longer being spherical. The fourth 
is the elastic strain energy, which depends quadratically on 
the difference between e, the actual shape of the star, and 
eo- The constant A depends on the stellar equation of state, 
and will be of the order of the gravitational binding energy 
of the star. The constant B also depends on the equation of 
state, and will be of order of the total electrostatic binding 
energy of the ionic crustal lattice. The equilibrium configura- 
tion can be found by minimising the energy at fixed angular 
momentum: 

0€ J 

For an entirely fluid star we would put -5 = 0, giving an 
oblateness of order of the ratio of kinetic and gravitational 
energies per unit mass: 



Efluid 



T 

-'star ^ ^ 

AA ■ 



(3) 



Given that A is of the order of the gravitational binding 
energy we can write this as: 



Efluid 



GM 



(lok) 



(4) 



where ft = 2nf is the angular frequency, -Re the neutron 
star radius in units of lO^'cm, and Mia the mass in units of 
1.4M0. 

When the strain term is included we find 



/starfi^ B 
e — -rr-. + . , .-, £0 = en -f oeo. 



4(A + B) A-^B 



(5) 



The oblateness is made up of two parts. The first, en, scales 
as and is due to centrifugal forces. We will refer to this 
as the centrifugal deformation. The second term, feeo, is due 



entirely to the stresses of the crystalline solid, and will be 
referred to as the Coulomb deformation. We have defined b = 
B/{A + B), which we will refer to as the rigidity parameter. 
It is equal to zero for a fluid star (-B — 0) and unity for 
a perfectly rigid one {B/A — + oo). Realistic neutron star 
equations of state imply that b takes a value of: 



1.6 X IQ-^Rl/Mi 



(6) 



(See Jones (2000) for a simple derivation, and Ushomirsky, 
Cutler & Bildsten (2000) for a detailed numerical treat- 
ment). Because this is small, b is approximately equal to 
B/A, and so is simply the ratio of the crustal electrostatic 
binding energy to the total stellar gravitational binding en- 
ergy. It is this second deformation that makes free precession 
possible, not the (possibly much larger) centrifugal deforma- 
tion. We will write the total change in the moment of inertia 
tensor due to this Coulomb term as A7d, so that 

A7d 



= feeo. 



(7) 



Thus we see that the effect of elastic stresses in the 
crust is to change the shape only slightly from that of the 
corresponding fluid body. Physically, the smallness of this 
distortion is due to the Coulomb forces being much smaller 
than the gravitational and centrifugal ones. 

As b is small we have e « en ~ efluid, so that in the 
following sections we will use equation (^) to estimate a 
star's actual oblateness. Also, e and eo can differ at most by 
the breaking strain Ubroak of the crust, so we expect: 



|e — eo| < Ubrcak. 



(8) 



This breaking strain is very poorly constrained. Estimates 
have ranged from ~ 10~^ to values very much lower — see 
section ^. 

In the following section it will be useful to have an ex- 
pression for the ratio of the crustal to total stellar moments 
of inertia. From Ravenhall & Pethick (1994) we flnd: 

-/crust 



Ista 



1.5 X lO'^Rs/Mt.i- 



(9) 



We will approximate the total stellar moment of inertia us- 
ing the constant density result 

/star = ImR" = 1.12 X lO^"" gem'' M^aRI (10) 
5 



3 MODELLING FREE PRECESSION 

Real neutron stars will consist of an inelastic crust contain- 
ing, in a non-spherical cavity, a compressible liquid core. 
This core will be made up of a viscous electron-proton 
plasma, a neutron su perfluid, and poss ibly more exotic 
phases of matter also (Glendenning 1997). The crust itself 
may contain a pinned superfluid in its inner parts. Also, a 
magnetic field will thread the star, in a way which depends 
on the properties of the superfluid phase. The superfluid 
core will couple in a frictional way to the crust. The free 
precession of such a system will clearly be far more com- 
plicated than that of a rigid body. The strategy that has 
been employed to explore this free precession is to look at 
only one complicating factor at a time. Following this ap- 
proach, we will briefly describe each complicating feature. 
In this section the effects of elasticity, crust-core coupling 
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then gives 




Figure 1. This figure shows the reference plane for a freely pro- 
cessing body, which contains the deformation axis n^j , the angular 
velocity vector f2 and the fixed angular momentum J. The vectors 
and f2 rotate around J at the inertial precession frequency <f>- 
We refer to 9 as the wobble angle. 



and superfluid pinning are described, and then combined 
to give a reahstic free precession model. All of these com- 
plications preserve the kinetic energy of a precessing body. 
The dissipative effects of a frictional crust-core coupling and 
gravitational radiation reaction are described in section U. 



3.1 Rigid shell 

We will begin by describing the simple case of a rigid shell. 
The moment of inertia tensor of any axisymmetric rigid 
body can be written as 



I = /o<5-f A/d(ndnd- V3), 



(11) 



where the unit vector nd points along the body's symmetry 
axis. Then the principal moments are Ii — I2 = lo — AJd/3, 
I3 = lo + 2A7d/3, so that h - h = A7d. The angular 
momentum is then related to the angular velocity according 
to 



J = (Jo - A/d/3)n - A/d^snd, 



(12) 



where the 3-axis lies along nd- This shows that the three 
vectors J, fl and nd are always coplanar. Following Pines & 
Shaham (1972) we will call the plane so defined the reference 
plane (see figure |l]). Given that the angular momentum is 
fixed, this plane must revolve around J. The free precession 
is conveniently parameterised by the angle 6 between nd 
and J. We will refer to this as the wobble angle. For a nearly 
spherical body the angle 9 between ft and J is much smaller 
than the angle between J and nd, according to 

e ^ sine* cos 6». (13) 
Ii 

We will denote by nj the unit vector along J. Decomposing 
the angular velocity according to 



A/d 



(15) 
(16) 



The symmetry axis nd performs s a rotation about J in 
a cone of half-angle 9 at the angular frequency 0. We will 
refer to this as the inertial precession frequency. There is 
a superimposed rotation about the symmetry axis nd at 
the angular velocity tp . This is usually referred to as the 
body frame precessional frequency, with the corresponding 
periodicity known as the free precession period: 

P.. - ^ (17) 

For a nearly spherical body equation (^6|) shows that ip <^ <j>, 
or equivalently P <C Pip. Note that the angles {9,(j),tp) axe 
simply the usual Euler angles which describe the orientation 
of the rigid body (see e.g. Landau & Lifshitz 1976, figure 47). 

In subsequent sections we will find it necessary to make 
the approximations of small wobble angle and nearly spher- 
ical stars. In this case equations ([!) and (|l|) become: 

0^^9, 
Jo 



(18) 



A/d 



(19) 



3.2 Elastic shell 



Following section g, we will write the moment of inertia ten- 
sor of a rotating elastic shell as the sum of a spherical and 
two quadrupolar parts: 



1 = I0S + A/d(ndnd - S/3) + A/n(nnno - S/3). 



(20) 



(See Pines & Shaham (1972) and Munk & McDonald (1960) 
for further discussion) . The first term on the right hand side 
is the moment of inertia of the non-rotating undeformed 
spherical shell. The second term is the change due to crustal 
Coulomb forces, and has the unit vector nd, fixed in the 
crust, as its symmetry axis. The third term is the change 
due to centrifugal forces, and has nn, the unit vector along 
ft, as its symmetry axis. When ft moves with respect to 
the body the shell changes shape. This is why the shell is 
described as elastic. 

When the directions nd and no coincide the body spins 
about its symmetry axis without precessing. When the two 
directions differ the body will precess with triaxial shape. 
Proceeding exactly as in the rigid body case, it is easy to 
prove that the elastic body undergoes a free precessional 
motion like that of a rigid axisymmetric body, with what 
we will call an effective moment of inertia tensor given by: 



h = Io~ AId/3 + 2A/n/3, 
73 = /o + 2A/d/3 + 2A/n/3. 



(21) 
(22) 
(23) 



ft 



(14) 



Crucially, the centrifugal piece of the quadrupole moment 
enters in a spherical way — it is only the A/d piece that is 
responsible for the free precession. In particular, equations 
( p^ and (^) still apply, with /i and I3 as given above, and 
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with the precession frequency and angle 6 both proportional 
to the Coulomb deformation A/d. 



3.3 Rigid shell containing non-spherical fluid 

First consider the free precession of a rigid axisymmetric 
shell containing a spherical fluid cavity. We will assume that 
there are no viscous or frictional interactions between the 
shell and the fluid. In this case the free precession of the shell 
is exactly the same as if the fluid were not there. However, 
if the c avity is non-spherical, there will be a reaction force 



between the rigid shell and the fluid, due to the fluid tend- 
ing toward a configuration symmetric about its rotational 
axis, and therefore 'pushing' on the shell. This pushing is 
known as inertial coupling. In the case of a homogeneous 
incompressible fluid, the combined motion of fluid and shell 
is given in Lamb's monograph (1952). A number of simpli- 
fying assumptions are necessary: The motion of the fluid is 
always one of uniform vorticity; the ellipticity of the shell 
and cavity are small; and the wobble angle 6 is small. 

It is then found that the shell undergoes the usual free 
precession motion, so that equations ( p^ and ( p^ ) apply, 
with A/d equal to the difference between the 1 and 3 princi- 
pal moments of inertia of the whole star, not just the shell, 
while Jo refers to the shell only. The system's total angular 
momentum J is simply the sum of the angular momenta of 
the shell and fluid, both of which remain very nearly parallel 
to J. As defined previously, 9 is the angle that the angular 
velocity of the shell makes with J. Crucially, the angular 
velocity of the fluid remains very close to the fixed total an- 
gular momentum of the system, i.e. inertial coupling does 
not cause the fluid to participate in the free precessional mo- 
tion of the shell. 

3.4 Rigid shell with pinned superfluid 

According to many neutron star models which attempt to 
explain post-glitch behaviour, the neutron vortices which 
coexist with the inner crust become pinned, at many points 
along their length, to nuclei in the crustal lattice. The veloc- 
ity field of this pinned superfiuid is specified entirely by the 
distribution of these pinning sites. The pinning sites them- 
selves are rigidly fixed to the crust. It follows that if a rigid 
crust is set into free precession, the instantaneous velocity 
field of the pinned superfluid continually adjusts, according 
to the orientation of the crust. Such a precessing star was 
studied by Shaham (1977). If we write down the angular mo- 
mentum of the crust, including that of its pinned superfiuid, 
we have 



J = in + Js 



(24) 



The quantities I and f2 refer to the crust only, while Jsf 
is the angular momentum of the pinned superfiuid. The ori- 
entation of Jsf is fixed with respect to the crust. We will 
consider only the simplest case, where all the pinned su- 
perfluid points along the crust's deformation axis, s o th at 
Jsf ~ Jspnd. Then repeating the analysis of section 3.1 we 
find (for nearly spherical bodies with 6^ ^ 1) 

A/. Jsf 



(25) 



AJ 
lo" 



Jsf 



(26) 



In words, the pinned superfluid acts so as to increase the ef- 
fective non-sphericity of an oblate body, increasing the mag- 
nitude of both the body frame precession frequency, and the 
misalignment between the crust's angular momentum and 
angular velocity vectors. 

The rotation rate of the pinned superfluid will prob- 
ably be very close to the star's rotation rate, as even a 
small difference between the two wou ld give rise to a Magnus 



force which would break the pinning ( Lyne & Graham- Smith 
199S). We can therefore put Jsf ~ /sf0, where Jsf is the 
moment of inertia of the pinned superfiuid. Then the above 
equations can be written very simply as: 



Aid 
Jo 

Aid Jsf 
Jo Jo 



Jsf 



(27) 



(28) 



We will call the term in square brackets the effective oblate- 
ness parameter. It is made up of both crustal distortion and 
pinned superfluid parts. 

3.5 Composite model 

We now wish to consider the free precession of a more realis- 
tic composite model — an elastic shell containing a fluid cav- 
ity, with superfluid pinned to the shell. We would then form 
the equation of motion of the body by combining equations 
(|^ and When inertial coupling forces are neglected it 
is straightforward to repeat the analysis to show that the ef- 
fects of elastic deformation and superfluidity add in a simple 
way. (Equations (|2^)"(|2^) apply, with AJd the Coulomb de- 
formation of the shell and Jq the crustal moment of inertia) . 
The effect of also including inertial coupling forces will be 
to set the AJd factor to the deformation in the moment of 
inertia of the whole star, not just the change in the crustal 
moment of inertia. Then equations of the form ([27[)-(|2^) 
apply again. In full: 



= ecftO, 
AJd 



(29) 
(30) 

(31) 



and 



Jsf 

where Jo is equal to the crustal moment of inertia only, while 
AJd is the Coulomb-induced deformation in the moment of 
inertia of the whole star, and Jsf is the moment of inertia 
of the pinned superfluid. 



4 DISSIPATION MECHANISMS 

The model described above would, once excited, precess for- 
ever, as no dissipative energy losses have been included. A 
real star will suffer a number of such losses. To complete our 
model we will therefore consider two types of dissipation: A 
frictional crust-core coupling, and gravitational radiation re- 
action. 

Before considering these particular cases we will derive a 
general expression for the damping timescale. The energy of 
the precessing crust can always be written as a function of its 
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angular momentum and wobble angle. For wobble damping 
the angular momentum of the crust is nearly constant, so 
we can write: 
idE 



89 



(32) 



Here E denotes the total energy of the shell plus pinned 
superfluid, and will include kinetic, elastic and gravitational 
parts. Cutler & Jones (2000) have shown that, to leading 
order, only the kinetic energy need be considered. It is then 
stra ightforward to take the free precession model of section 
and evaluate the partial derivative of the denominator. 



3.5 



givmg 



E 



Expressed as a damping time this is 



re 



(33) 



(34) 



The quantity E will always be negative, corresponding 
to the conversion of mechanical energy into heat or radia- 
tion. It therefore follows that dissipation tends to decrease 
the wobble angle of stars with oblate deformations, but tends 
to inc rease the wobble angl e of stars with prolate deforma- 
tions (Cutler & Jones 2000). The crustal strains considered 



in this paper will almost certainly lead only to oblate de- 
formations. More exotic scenarios (perhaps a very strong 
toroidal magnetic field) might lead to prolate ones. Dissipa- 
tion in such stars would eventually lead to the deformation 
axis being orthogonal to the spin axis. Such a non- precessing 
triaxial star would then spin down gravitationally. 



4.1 Frictional crust-core coupling 

We have already considered one crust-core interaction, 
namely inertial coupling, due to the core fluid simply 'push- 
ing' on the precessing shell. However, in a real star there will 
be additional crust-core interactions, which have been inves- 
tigated by theorists attempting to explain the post-glitch 
behaviour of pulsars. The core itself consists of two coexist- 
ing fluids: a plasma of electrons and protons, and a neutron 
superfluid. The plasma makes up only a few percent of the 
total mass. 

The crust-core interaction proceeds in two stages. In the 
first stage the crust couples to the core plasma. This cou- 
pling is mediated by two separate interactions: plasma vis- 
cosity, and an electromagnetic coupling, where Alfven waves 
communicate the crust's motion to the interior. In the sec- 
ond stage, the plasma couples to the neutron superfiuid, 
due to the scattering of electrons off the superfiuid vortices. 
(Strictly, the electrons are scattered by the magnetic field 
created by protons which are entrained around the vortex 
cores; see Alpar & Sauls 1988). This second interaction is 
frictional in nature, i.e. is a drag force proportional to the 
velocity difference between the electrons and the vortices. It 
is sometimes referred to as 'mutual friction'. 

There will be energy losses associated with both cou- 
plings. However, following the work of Easson (1979) it is 
usually assumed that the crust-plasma coupling timescale 
is much less than the frictional one, so that the crust and 



core plasma can be treated as a single system, interact- 
ing frictionally with the neutron core. We will work in this 
limit. (Relaxation of this assumption will in fact tend to 
strengthen our conclusion, namely that the neutron super- 
fiuid does not follow the precession of the crust). 

This frictional interaction is important in two regards. 
Firstly, it will lead to a dissipation of the precessional energy. 
Secondly, if strong enough, it would cause the whole star — 
crust and core — to precess as a single body. We can therefore 
imagine two extreme cases. When the frictional coupling is 
very weak the crust precesses on top of a no n-precessing core, 
with only the inertial forces of section 3.3 coupling the two. 
Then equation (^l|) applies, with 7o = /crust. At the other 
extreme, when the frictional forces are very strong, the star 
precesses as a single unit. Then equation (|3l| ) again applies, 
but with Jo ~ Jstar. It is clearly important to know where 
on this scale real neutron stars can be expected to be found, 
both from the point of view of modelling the precession, and 
for estimating the rate at which it is damped. 

The motion of a rigid shell containing a spherical cav- 
ity, with a frictional coupling acting between the two, was 
invested by Bondi & Gold (1955). They made the simplify- 
ing assumption that the motion of the fiuid was one of rigid 
rotation. The shell is acted upon by a torque 



T = ^^^(fjfluid — nsolid) 



(35) 



where K is & positive constant. An equal and opposite torque 
acts on the fluid. The equations of motion are then 



dt 



A(^ifluid — S^soUd) 



dt 



(36) 



Bondi & Gold found the normal modes of this two compo- 
nent system for small wobble angles. The results below can 
be readily found from their solution. 

First consider the simple case of the non-precessing shell 
and fluid, rotating about the same axis, but at different 
rates. Using the component of the above equation along the 
shell's symmetry axis, it is easy to show that the relative 
rotation rate decreases exponentially. We will write the e- 
folding time of this decay as n rotation periods, correspond- 
ing to a time 2Tin/Q.. 

Using the components of ( ^ij ) orthogonal to the sym- 
metry axis it can then be shown that in the case n the 
fluid is tightly coupled to the shell, and such a body would 
rotate and precess as if it were a single uniform solid. In the 
case n 2> 1 the fluid is only loosely coupled, so that while the 
shell undergoes free precession, the fluid's angular velocity 
vector remains very nearly flxed in space. 

Alpar, Langar & Sauls (1984) calculated the coupling 
strength for the electron-vortex core interaction described 
above. Using this result, Alpar & Sauls (1988) estimated n 
to lie in the interval 400 10^. As n 1 it follows at once 
that real neutron stars lie in the weak coupling regime. In 
this regime it can then be shown (using the components of 
( p6| ) orthogonal to the symmetry axis) that if the shell is 
set into free precession, with its angular momentum vector 
remaining along the spin axis of the fluid, the free preces- 
sion is damped, with an e-folding time of n body frame free 
precession periods, i.e. in a time 27rn/i/). 

In summary, this frictional interaction is too weak to 
cause the core neutron superfluid to participate in the free 
precession of the star's crust. Instead, it serves only to damp 
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the free precession of the crust, on a timescale of between 400 
and 10* free precession periods. 



4.2 Gravitational radiation reaction 

The effect of gravitational radiation reaction on precessing 
elastic bodies with a fluid non-spherical core was recently 
considered by Cutler & Jones (2000). It was found that the 
wobble angle was damped exponentially on a timescale 



5c' 



1 



2G (A7d)2 n*' 



(37) 



This remains valid when the effects of superfluid pinning are 
included. Parameterising: 



1.8 X 10 yr 



10''^gcm2\ / 10" 



VA/d//o 



/ \ 1ms/ 



For instance, for a star with A/d//o ~ fcff = 10~^, lo ~ 
10'*^gcm'^ and P — 1ms, this corresponds to damping in 
n = 6 X 10** free precession periods. Comparing this with the 
damping rate due to friction, we see that gravitational radi- 
ation reaction is not an important source of wobble damping 
in any neutron star of physical interest. (This does not mean 
that the spin-down component of the radiation reaction is 
unimportant — see section [7.2| ). 



5 CRUST FRACTURE AND MAXIMUM 
WOBBLE ANGLE 

A real neutron star crust will have a finite breaking strain 
Wbroak • The actual value of this breaking strain is highly un- 
certain. As discussed by Ruderman (1992), extrapolations 
from laboratory crystals suggest that values as high as 10~^ 
are possible, but the actual value may be much lower. Rud- 
erman suggests a value of 10~* as plausible. For an axisym- 
metric star this breaking strain will limit the difference be- 
tween its actual shape, which is approximately efluid, and the 
reference shape of the crust, eo- However, when a star is set 
into free precession, additional time-dependent strains will 
result, even if the star is initially 'relaxed' (i.e. efluid = eo). 
In gene ral, a precessing star's crust will suffer both sorts of 
strain (Pines fc Shaham 1972). 



First consider a non-precessing spirming-down star. Its 
crust presumably solidified (and eo was fixed) when it was 
spinning more rapidly that at present. This means its refer- 
ence oblateness eo is greater than its current actual oblate- 
ness efluid, so that we would now expect eo > efluid. Combin- 
ing with the bound due to crust fracture we then have 



< eo — efluid < Itbroak- 



(39) 



Even if the crust has cracked a number of times during spin- 
down we would still expect the star to have some 'memory' 
of its initial shape (unless the cracking was able to relieve all 
the stresses in the crust), so that the above equation should 
still hold. 

In the case of an accreting star the situation is different. 
If the star has been spun up from a relaxed state at a small 
rotation rate we would expect its current oblateness efluid to 
exceed its reference oblateness eo. However, the accretion will 
tend to create new crustal material relaxed to the current 



rotation rate of the star, so that the difference efluid — eo may 
be rather small. 

We will now look at the opposite case, when the strains 
are due to precession only, i.e. when eo = efluid. A simple 
treatment is possible, based upon the known geometry of 



free precession for an elastic body, as discussed in section 3.2 



These strains can be used to place a limit on the maximum 
possible wobble angle 9 a star can sustain. 

We know that (for small wobble angles at least) when an 
initially axially symmetric body is set into free precession a 
deformation Aid remains along the axis nd fixed in the star, 
while a deformation A/n points along the angular velocity 
vector. From the point of view of an observer attached to 
the crust, a deformation of size AIn describes a cone of 
half-angle 9 + 6^6 about nj. This change in shape is 
all we need to know to estimate the strain: The change in 
position of any given particle is of order Ren9, while the 
corresponding strain is of order enS. This precession-induced 
strain is not constant, but varies with magnitude en9 over 
one (body frame) free precession period. As there exists a 
maximum strain Ubroak that the solid can withstand prior 
to fracture, the wobble angle will be limited to a value of 

^tbroak/en ~ "Ubrcak/efluid SO that: 



0.45 



100 Hz 



^brcak 

10-3 



radians. 



(40) 



Qualitatively, we can say that faster spinning neutron stars 
have larger bulges to re-orientate and therefore can sus- 
tain smaller wobble angles prior to fracture. For sufflciently 
slowly spinning stars the above equation breaks down, yield- 
ing angles in excess of 7r/2. The wobble angles of such slowly 
spinning stars are not limited by crustal strain. To give two 
extreme examples, for a breaking strain of 10~^ the wobble 
angle of a 300Hz neutron star in a low-mass X-ray binary 
or millisecond pulsar would be limited to about 3°, while a 
'standard' field pulsar spinning at around a Hz could precess 
with any wobble angle. 

We therefore have two separate bounds that a freely 
precessing star must satisfy, one due to the mismatch be- 
tween its reference and actual shape, and one due to free 
precession. These bounds will be of use in section ^, when 
we examine possible observations of free precession. 



6 EFFECT OF FREE PRECESSION ON THE 
PULSAR SIGNAL 

In this section we will examine the effect of free precession on 
the electromagnetic signal of a pulsar. The aim is to describe 
how free precession might be detected using electromagnetic 
data, by examining variations in the pulse frequency, am- 
plitude and polarisation. From these we hope to test our 
free precession model and extract information concerning 
the wobble angle, effective oblateness and superfiuid pin- 
ning. 

we will consider the problem of how the 



In section 3.1 



precession affects the electromagnetic pulses for a torque- 
free top. This calculation requires only the geometry of free 
precession, and we will refer to these modulations as the 
purely geometric modulations. However, real pulsars are 
acted on by electromagnetic spin-down torques. The magni- 
tude of these torques will themselves be modulated by the 
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free precession, although the form of the modulation de- 
pends on the model of spin-down torque employed. T hese 



torques will, in turn, modify the pulsations. In section 
we will include these variable torques using a simple spin- 
down model. 

Some of the results given here have been presented pre- 
viously, in varying degrees of generality. We hope that by 
assembling all of the important observational characteristics 
of free precession in one place, we will provide a resource for 
observers wishing to asses the likelihood of a given pulse 
modulation being caused by free precession. 



6.1 Effect of free precession on the pulsar signal: 
Geometric effects 

We will begin by modelling the neutron star as a torque-free 
precessing symmetric top. We will model the pulsations in 
the obvious way, i.e. as a thin beam fixed with respect to the 
star's body axes, aligned with the dipole moment m. Each 
pulsation then corresponds to the passage of m through the 
plane containing the angular momentum vector J and the 
observer. The motion of m is then a slow rotation at ip 
about nd, the deformation axis, superimposed on the rapid 
rotation of nd at <f> about J. This will lead to variations in 
the pulse phase, amplitude and polarisation, on the body 
frame precession timescale Pfp = 2n/^. We will begin by 
considering the phase variations. 



6.1.1 Phase modulation 

The problem of phase modulation due to precession was 
first considered by Ruderman (1970), and has been elab- 
orated upon since by Bisnovatyi-Kogan et al. (1990) and 
Bisnovatyi-Kogan & Kahabka (1993) in connection with the 
35 day periodicity observed in Her X-1. Following the dis- 
covery of planets around pulsar PSR 1257-1-12 a number of 
authors re-examined the issue to check whether or not free 
precessi on could mimic planetary perturl a ations of the pul 



sation (Nelson, Finn fc Wasserman 199o| ; [Cordes 1993 



et al. 1993 



Gil 



findings in section 3.2 



Glendenning 199£| ) . We will comment upon their 



Let m denote a unit vector along m. Denote the ori- 
entation of the body frame axes {x, y, z} with respect to 
the inertial frame axes {x, y, z} by the usual Euler angles 
{6,<j>,'<P). (See figure 47 of Landau & Lifshitz 1976). Let ih 
lie in the xz plane, at an angle x < ''r/2 to the i-axis. Then 
the components [mx,my,mx] of m, referred to the inertial 
frame are 

cos (j) cos tl) sin x — sin <j) cos 6 sin sin x + sin sin 9 cos x 
sin (j) cos if) sin x + cos (p cos 6 sin tp sin x — cos <p sin 9 cos x 
sin 9 sin ^p sin x + cos 9 cos x 

Define Euler-like angles Q and $ to describe the polar angle 
and azimuth of m. Then the $ angle describes the phase of 
the pulsar signal, with a pulsation being observed whenever 
$ is equal to the azimuth of the observer, e.g. whenever 
$ = for an observer in the inertial a; > 0, 2 > quarter- 
plane. Then 



tan $ = —2- . 



(41) 



and also 



+ arctan 



1 



cos tp tan X 



cos 9 \ tan 9 — sin ip tan x 



(42) 



(43) 



+tp sin X 



cos 9 sin X — sin tp sin 9 cos x 



(sin 9 cos X — cos 9 sin -ip sin x)^ + cos^ tp sin x 



where $ is the instantaneous electromagnetic frequency. Its 
time-averaged value is what we would normally refer to as 
the 'spin frequency' of the star, which we will denote by fi. 
In order to proceed it is necessary to treat the 9 > x a-nd 
9 < X cases separately. 



9 > X case 

First consider the average electromagnetic pulse frequency. 
An increase of (p by 27r at fixed tp causes m to rotate once 
about J, i.e. causes a pulsation. However, an increase in tp by 
27r at fixed 4> does not rotate m about J, i.e. does not cause 
a pulsation. It follows that the average electromagnetic fre- 
quency is exactly <p. The departure from this average spin 
rate can best be expressed as a phase residual A"l>: 



A* = $ 



arctan 



1 



cos tp tan X 
cos 9 \ tan 9 — sin tp tan x 



(44) 
(45) 



The denominator of the term in curved brackets is clearly 
non-zero for all tp, and so remains in the r ange —n/2 < 
A<I> < Tv/2 ( Nfelson, Finn fc Wasserman 199C| ). In the case 
S S> X wc find 



A<E> = cos tp 
sm9 



and 

A<i> : 



-tP- 



■ sin tp 



(46) 



(47) 



When X = 0, as would be the case for a star whose deforma- 
tion is due entirely to axisymmetric magnetic stresses, there 
is no modulation in the pulsations at all. The free precession 
of such a star would only be detectable is there was some 
non-axisymmetry in the pulsar beam. 



9 < X case 

As described in section ^ the wobble angles of rapidly ro- 
tating stars are limited to small values by the finite crustal 
breaking strain, so that for such stars this is almost certainly 
the case of interest. Again begin by considering the average 
pulsation frequency. An increase of (p by 2n at fixed tp causes 
rh to rotate once about J, i.e. causes a pulsation. It is also 
the case that an increase of tp by 27r at fixed (p causes m to 
rotate once about J, i.e. causes a pulsation. It follows that 
the average pulsation frequency is (p + tp. The phase residual 
is now given by 



A little algebra leads to 



A* = $ 



arctan 



-(<^ + ^) (48) 

(cos ^ — 1) sin tp sin x — sin 9 cos x 
cos tp sin X + (cos 9 cos tp sin x — sin 9 cos x) tan tp 
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It is straightforward to show that the denominator of this 
function is non-zero for all and so A$ remains in the range 
-7r/2 < A$ < tt/2 ([Nelson, Finn fc Wasserman 1990|). 



When S <C X find 



and 



^cosx 
sinx 



cos tp 



(49) 



(50) 



. i ■ cos Y . 

A$ = ?/)6l smtp. 

smx 

We therefore see that for small wobble angles the phase 
residual varies sinusoidally on the (long) free precession 
timescale, with an amplitude 9. The fractional variation in 
pulsation frequency is of order Otp/fp ~ 9ecS- 



6.1.2 Amplitude modulation 

The precessional motion will modulate the amplitude of the 
pulsar signal, although the precise form of the modulation 
will depend upon the geometry of the emission region. As- 
suming emission axisymmetry about the dipole axis and that 
the intensity of emission falls off over an angular width W , 
the fractional change in amplitude l\A /A due to precessional 
modulation is of order (Cordes 1993) 



AA 



Ae 

w ' 



(51) 



where AG denotes the change in polar angle of m over a 
free precession period. We have 



cos = = sin 6 sin ip sin x + cos 6 cos x, 



(52) 



from which we find the (obvious) results that Q has a max- 
imum value of X + ^ and a minimum value of |x — It 
follows that A0 is approximately equal to 6 or x, whichever 
is smaller. 

Pulsar duty cycles (i. e. angular widths of the beam) are 
typically of order of 10° (Lyne & Graham-Smith 1998), so 



we will parameterise according to 



M^exio- 

A 



(lO-O 



10° 

W 



(53) 



assuming 9 < x- This modulation will occur at the body 
frame free precession frequency. 



6.1.3 Polarisation modulation 

As above, we will assume a pulsar beam structure symmet- 
ric about m. In addition, we will take the polarisation model 
described in Lyne & Graham-Smith (1998, section 12.2). In 
this model the polarisation vector of the observed radiation 
is parallel to the magnetic field line where the radiation was 
produced. Then the time variation of the linear polarisation 
angle. A, can be calculated (see their equation 12.2). It is 
this quantity that is measured by observers. In the absence 
of free precession this angle varies rapidly at the spin fre- 
quency. Free precession will cause an additional modulation 
at precession period Pfp. A useful diagnostic for observers is 
the maximum rate of change of this angle with time, Amax, 
which occurs when the dipole axis lies closest to the ob- 
servers line-of-sight. Let i denote the inclination angle, i.e. 



the angle between J and the line-of-sight. Using the equa- 
tion of Lyne & Graham-Smith for A we then find (for 9 < x) 
a fractional modulation in Amax of 



AA„ 



2 sin i 



sinxsin(i — x) 



(54) 



6.2 Effect of free precession on the pulsar signal: 
Electromagnetic torque effects 

The calculations of section considered a precessing sym- 
metric top entirely free of torques. However, the fact that 
the star is visible as a pulsar means that it will be acted 
upon by an electromagnetic torque, and this should be in- 
cluded in the calculation. This torque will not affect the 
amplitude and polarisation arguments. However, provided 
that the torque is a function of the spin rate and orientation 
of m, i.e. of <& and O, the spin-down torque will be modu- 
lated by the precession. This modulation in the torque must 
be taken into account when calculating the phase residu- 
als. We will call this electromagnetic modulation. The affect 
on the phase residuals of this varying torque has been con- 
sidered analytically by Jones (1988) and Cordes (1993) for 
general torque functions, and numerically by Melatos (1999) 
for the Deutsch (1955) torque. We, however, will give a sim- 
ple description based on the vacuum point-dipole spin-down 
torque so that 



<E> = sin^ e. 



(55) 



where A; is a negative constant. The fractional change in 
spin-down rate due to precession is given by 



A£ 
IT 



A$ „A(sine) 



sin O 



(56) 



The prefix A denotes the departure of the respective quan- 
tity from the smooth non-precessing power law spin-down. 
We will consider the case 9 < x- From equation ( pOj) w e see 
that the first term is of order ecs9. Using equation (|52|) it is 
easy to show that 



sin 



sin^ X - 



29 sin X cos x sin t/j 



(57) 



so that the second term of ( |56[ ) is of order 9, and therefore 
is the dominant one. We then have 



A$ ~ — 2fcr2 sin X cos X sin 



(58) 



All the quantities on the right hand side are constant, apart 
from the angle ip, which is a linear function of time: ip = tpt. 
We can integrate once to get the perturbation in frequency 



A<i> ~ 2fcf2 sinxcosxcos?/> — , 

and once more to obtain the phase residual 

3 d 
A<[> ~ 2fcf2 sinxcosxsinf/;-:— . 

Now use iV'/'i'l = Ecff to give 

, ^ 1 d P 
A* = -cotx— , 



(59) 



(60) 



(61) 



where P denotes the spin period and Te = \^/^\ is the 
spin-down timescale. For example, if we consider a star 
where superfiuid pinning is not operative, and where only 
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the crust participates in the free precession, we can put 

Ecff ~ feen /star //crust (aSSUHliug 60 ~ Efluid) tO givC 



A$ = 25cotx(^) 



10" 



30 ms 



(62) 



lO'^ yrs 



t//s, 



0.015 



We have parameterised with a young pulsar in mind. This 
should be compared with the torque-free residual of equa- 
tion @. We see that for the parameterisation above the 
electromagnetic torque variation has greatly amplified the 
residual. Note, however, that for millisecond pulsars which 
have Tc 10^ years the torque variations are unimportant. 

Note also that for the parameterisations above, as A$ 
is greater than 2tt, this phase variation makes it difficult to 
detect such pulsars — any search algorithm that integrates 
radio data assuming a constant frequency source will go out 
of phase with the signal in an interval of the order of the free 
precession period. However, pulsar physicists do take such 
a modulation into account when performing data analysis 
when they search for binary pulsars, as a sinusoidal phase 
variation is also produced by a binary orbit: Small-angle free 
precession and nearly-circular binary orbits both produce si- 
nusoidal phase residuals, to leading order in wobble angle 
and orbital eccentricity, respectively. As pointed out by Nel- 
son et al. (1990), in order to distinguish between the two 
models it is necessary to include higher order terms. Then 
the residuals have different forms, allowing differentiation 
between the two models. Also, free precession will almost 
certainly lead to an amplitude modulation in phase with 
the timing residuals, whereas a binary companion would 
only produce amplitude modulation if it happened to cut 
the observer's line-of-sight onto the pulsar beam, providing 
another means of differentiation. 

We will present the phase residual in one more form. If 
we put ip — 27r/Pfp we obtain 



(63) 



7 ANALYSIS OF OBSERVATIONS OF FREE 
PRECESSION 

Having set out our free precession model in some detail, and 
described how free precession modulates the electromagnetic 
signal of a pulsar, we will now turn to the problem of extract- 
ing useful information from the pulsar observations. First we 
will assemble the necessary equations. 



7.1 Formulae required to extract source 
parameters from observations 

The problem divides into two parts: Extracting the wobble 
angle 9, and extracting information concerning the structure 
of the star, such as its reference shape or the amount of 
pinned superffuid. 

The extraction of the wobble angle is relatively straight- 
forward. We need only invert the equations of the previous 
section. Comparing the geometric phase residual of equation 
(Ml) with the electromagnetic torque residual of equation 



(|63|), we see that the electromagnetic torque significantly 
amplifies the geometric residual when 



> 1. 



(64) 



When this equality is satisfied equation (^3|) applies, and the 
wobble angle can be extracted (up to a factor of tan x) from 
the observed values of spin period, free precession period 
and phase residual magnitude according to 



P 
/^ 



A$tanx. 



(65) 



Note that if the spin-down torque is a steeper function of 
Q than the vacuum dipole model predicts (equation the 
wobble angle as estimated by the above equation will be an 
overestimate. If the torque is not as steep as assumed, then 
the calculated wobble angle will be an underestimate. 

When the inequality is reversed the phase residuals are 
described accurately by the geometric variation, and equa- 
tion (^) applies. Then the wobble angle can be extracted 
(up to a factor of tan x) from the value of A"l> according to 



= A$tanx. 



(66) 



Information concerning the structure of the star is more 
difficult to obtain. We hope to extract this information from 
equation (jsij), where the quantity ecff is simply the ratio of 
spin and modulation periods, P/Pfp. Ideally, we would like 
to extract the three quantities Aid (the deformation in the 
moment of inertia tensor caused by Coulomb forces), Isf 
(the moment of inertia of the pinned superfiuid), and lo 
(the part of the moment of inertia which participates in the 
free precession). Clearly, the problem is underdetermined. 
However, the quantity A/d//star is, for a given oblateness eo, 
constrained by the equation of state according to equation 
(^), so that: 



P , /star , /SF 

Eoff = = beo— + —. 



(67) 



The problem then becomes one of extracting eo, /sf, and 
/q. The quantity eo is, for a spinning-down star, almost cer- 
tainly greater than or equal to the actual oblateness efluid 
(given by equation ^ , and cannot exceed efluid by more than 
itbroak. AIso, if our Understanding of crust-core coupling is 
correct, /o is simply the crustal moment of inertia. Given 
these assumptions, the above equation can be used to place 
limits on /sf, the least certain of all the stellar parameters. 

However, we will employ a slightly simpler strategy, for 
the following reason: If, as some glitch theories require, a 
few percent of the stars' moment of inertia was made up of a 
superffuid pinned to the inner crust, the effective oblateness 
parameter would be of order unity, so that Pfp ~ P. If the 
observations below really do represent free precession, this 
prediction has failed completely. We will therefore proceed 
as follows. We will begin by setting /sf to zero, allowing us 
to extract a reference oblateness: 

eo^CcffT-y^ (/sf = limit). (68) 

/star 

In the case where /o = /crust we can use equations (^) and 
(|) to give: 



eo = 10 ecff— =- — 
Ha 



[ Isf ~ limit). 



(69) 
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This can then be tested against inequality (^9|). There are 
three possible cases. 

Ca se I: If the inequaUty is violated because eg < equid, 



our mo lei requires modification. The most likely modifica- 



tion is that more than just the crust participates in the 
motion, so that Jo > /crust- This would be a sign of stronger 
crust-core coupling than anticipated. Setting Isf to a non- 
zero value would only serve to increase /q. 

Case II: If the inequality is satisfied, then the observa- 
tion is consistent with superfiuid pinning not playing a role 
in free precession. The possibility that superfiuid pinning is 
making up some part of the total ees remains, but Isf/Io 
can be no larger than ~ tcff- 

Case III: The inequality is violated because eo — tfluid > 
Uhraak - Of coursc, librca k is uukuowu, but is surcly less than 



10 (Ruderman 1992). If the inequality is violated even for 



this large breaking strain, then superfluidity is playing the 
dominant role in determining the free precession frequency. 
In this case ecfi Isf/Io, or equivalently 



Isf 



-fu- 



( eo = limit). 



(70) 



Again parameterising with /q = /crust using equation (^) we 
obtain: 



= 1.5 X 10" ec 



ML 



;eo = limit). (71) 



T.2 Analysis of observations 

We will now apply our free precession model to the pro- 
posed observations of free precession to have appeared in 
the literature. The collection below does not represent all 
of the proposed candidates: We have not included observa- 
tions where only the (proposed) free precession timescale 
Pfp has been measured, and not the (average) spin period 
P. Also, we have not included sources where only one modu- 
lation cycle in the pulsations (or less) has been observed. In 
restricting our sample in this way, we will confine our atten- 
tion to sources where the evidence for a modulation in the 
pulsations is reasonably secure, and where the observational 
data is good enough for source parameters to be estimated. 



PSR B0531+21 (The Crab pulsar) 

Lyne, Pritchard and Smith (1988) observed a variation in 
the phase residual of the Crab pulsar of about 1.9 radi- 
ans, with period 20 months. Jones (1988) suggested free 
precession as the cause, and pointed out that the residuals 
would be dominated by the electromagnetic torque varia- 
tion. Equation (pSj ) then gives 6 ~ 5 x 10"^ tanx radians. 

Equation (pM gives a reference oblateness of 6 x 10~^, 
while equation (0) gives tfluid = 2 x 10"'*. This violates the 
inequality ( |39[ ) (case I) — the star would seem to have a ref- 
erence shape much less oblate than its current actual shape. 
Even if the whole star is assumed to participate in the free 
precession, so that Iq = /star, we find eo — 4.3 x 10~^ (equa- 
tion |6^), so that inequality ( p9| ) is still violated. Therefore, 
it is not possible to make this observation fit our free preces- 
sion model for any sensible stellar parameters. Given that 
the modulation has not been detected in subsequent obser- 
vations, it seems unlikely that free precession was detected. 



More recently, evidence has been presented for a 60s 
modulation in the Crab pulsa r's amplitude and phas e resid- 
ual, both in the optical band (^adez fc Galicic 1996a ; Sadez 
& Galicic 1996b|; |Cadez fc Galicic fc Calvani 1997] ). This is by 



far the shortest free precession period to have been proposed 
in the literature. Cadez & Galicic report a phase residual 
amplitude of A<E> = 1.2 x 10"^ radians. For Pfp = 60s elec- 
tromagnetic torque amplification is insignificant. Equation 
( p6[ ) gives 6 1.2 x 10^^ tanx radians. They also report a 
fractional amplitude modulation of 6 x 10^^. Equation ( |53| ) 
then gives « 1.1 x 10"*. These two 9 estimates would agree 
for X ~ 0.1 radians. 

In the absence of pinning, equation (^) gives a refer- 
ence oblateness of 0.55. This is a large value, violating ( |39| ) 
(case III) even for Ubrcak = 10~^. In other words, the refer- 
ence oblateness is too large to be accounted for by Coulomb 
deformation alone. It is therefore necessary to invoke super- 
fiuid pinning: Using equation ( |7l| ) , we see a fraction of order 
10~^ of the star's moment of inertia needs to be pinned to 
fit the data. However, a subsequent search (Golden et al 



200C) has failed to confirm the modulation, weakening the 
precession hypothesis. 



PSR B0833-45 (The Vela pulsar) 

Deshpande & McCuUoch (1996) have presented evidence for 
a 165d variation in the Vela's amplitude in the radio band. 
The fractional modulation is of order 1/2. They have also 
investigated the difference in time-of-arrival of the pulses at 
two different radio frequencies, and found a 330d variation. 
The most natural explanation of this latter variation, and 
probably therefore of the former too, would be connected 



with refractive scintillations (Lyne & Graham-Smith 199J) 
due to the inter-stellar medium. However, the authors sug- 
gested free precession as the cause. It is not easy to see how 
precession could cause the variation in the time-of-arrival 
difference between the radio frequencies. We will therefore 
concentrate on the amplitude variation. Equation (js^) gives 
^ ~ 0.1 radians ~ 6°. 

Assuming no superfiuid pinning, the 165d modulation 
in the Vela gives a reference oblateness of eo = 6 x 10~^, 
a factor of 4 less than the calculated actual oblateness, 
efluid = 3 X 10^^. We therefore see that the inequality ( |39[ ) 
is violated (case I) — the star is more nearly spherical than 
we would expect. In order to increase eo to its minimum ac- 
ceptable value of efluid, at least 6% of the star would need 
to participate in the free precession. Inclusion of superfiuid 
pinning increases this value. 



PSR B1642-03 

Cordes (1993) has reported evidence for 10^ d variations in 
the pulse shape and timing residuals of pulsar B1642-03, 
both in the radio band. Three cycles have been observed 
(see figure VI of his paper), although the profile is far from 
sinusoidal. The fractional pulse shape modulation was ap- 
proximately 0.05, so that equation (^) leads to 6 ~ 9x 10"'' 
radians. The phase residual amplitude is difficult to identify 
as it seemed to increase over the observation period, but a 
value A<I> « 0.025 is a reasonable average. This pulsar has 
residuals dominated by the electromagnetic torque variation 
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(equation |6^, so that equation (|65|) gives S ~ 8 x 10~^ tan x 
radians. These two estimates are consistent for x ~ 0.1 ra- 
dians. 

For the reference oblateness in the zero pinning hmit, 
equation (p9|) gives eo = 5 x 10"^, while equation M) gives 
Efluid = 10 , so that inequality (|3^) is satisfied (case II). 
Therefore, this observation is consistent with zero superfluid 
pinning. The calculated reference oblateness is very close 
to the actual shape, suggesting that the star is virtually 
unstrained. This could imply a very low crustal breaking 
strain of order 10~®. The pinned superfluid component can 
make up no more than 10~^° of the total stellar moment of 
inertia. 



PSR B1828-11 

Very recently a 1009d periodicity has been rep orted in the 



1 p2 

A$ = %-e'^cos2tb 



(75) 



radi o sh ape and residuals of pulsar B 1828- 11 ( stairs ct al _ 
200C |). j Vpproximately four cycles have been observed. The 
pulse width is about 3°, and the pulse shape changes dras- 
tically over one free precession period, suggesting 6* 3°. 

The phase residuals give, via equation (|65|), a wobble 
angle of only 3.5 x 10~* radians = 0.02°, two orders of 
magnitude less. However, as is immediately apparent, there 
is a very strong 504d periodicity in the data. If we put 
Pip — 504d we obtain a, 6 — 0.08°, still a factor of order 
30 too small to agree with the wobble angle derived from 
the amplitude modulation. (Stairs et al. obtain 0.3° using a 
similar method of calculation, still one order of magnitude 
smaller than the amplitude-derived value). There is clearly 
a problem in extracting the wobble angle of the source. 

However, the fact that the data contains a strong peri- 
odicity at 504d strongly suggests the following: The free pre- 
cession period is indeed 1009 days, but the magnetic dipole 
lies very nearly orthogonal to the star's deformation axis, 
i.e. X ~ 7i"/2. Then both the phase modulation of equation 
( ^sl ) and the amplitude modulation connected with equa- 
tion (|5^) have significant components at 2i{j- To show this 
expand these equations to second order in 9. For the phase 
modulation: 



A* 



' cot X cos ■ 



■ sin 2^(1 + 2cot x)> 



(72) 



where the time dependence of the right hand side is due to 
ip = ^t. The first term is the linear (in 6) phase modulation 
at tp, and the second the quadratic (in 6) phase modulation 
at 2tp. The amplitude modulation for an unmagnetised star 
is due to the variation in O given by 



sin^ = sin^ X + cos 2x 



2 sin X cos x sin i^O 



— sin^ X cos 2i/). (73) 



The first two terms are uninteresting constants, the third 
the linear phase modulation at t/i, and the last the quadratic 
phase modulation at 2tp- The quadratic terms in (^^ and 
( [73I ) dominate the linear ones when: 

4 



tanx > 



(74) 



The calculation of section for the phase modulation for a 
magnetised star can then be repeated. In the x ~* '""/S limit 
we obtain 



We can then invert this relation to give the wobble angle. 
Using Pfp gives 9 = 2° ,in excellent agreement with the value 
estimated using the arnplitude modulation. For consistency 
X must then satisfy ([74|), which gives x > 89°. In other 
words, for this scenario to apply, we require near perfect or- 
thogonality between the deformation axis and the magnetic 
dipole axis. 

The spin and modulation periods of this star are, by co- 
incidence, almost identical to those reported above for PSR 
B1642-03, so the same conclusions regarding the star's struc- 
ture can be drawn, viz that the observation is consistent 
with our free precession model, with zero superfluid pinning 
and a crustal strain of no more than ~ lO^*". The maximum 
amount of superfluid pinning allowed is of order 10"^^°. 



Remnant of SN 1987A 

Middleditch et al. (2000a,b) have recently presented evi- 
dence for a 2.14ms optical pulsar in the remnant of SN 
1987A. The source was observed intermittently between 
1992 and 1997. A modulation in the phase and amplitude 
were detected, which the authors suggested may be due to 
free precession. The period of the modulation seemed to 
vary, spanning a range of 935s to 1430s. The amplitude of 
the phase modulation seemed to vary from ~ 48° to ~ 60°, 
although in some observing runs it was possibly zero. Mid- 
dleditch et al. also suggest that the spin-down of the pulsar 
is due to the gravitational radiation reaction associated with 
the free precession. 

We will examine the free precession and gravitational 
wave driven spin-down using our model. Crucially, the grav- 
itational wave spin-down hypothesis can be tested: The ratio 
of the spin to modulation periods will enable us to calculate 
the deformation Aid of the star, while the phase modula- 
tion will allow us to calculate the wobble angle 9. These can 
then be combined to give the gravitational wave spin-down, 
which can then be compared to the observed value of the 
frequency derivative. 

First we will extract a value for the wobble angle using 
the amplitude of the phase residual, A$. This amplitude is 
highly variable, implying a variable wobble angle. At times 
the phase modulation seemed to disappear. If real, this dis- 
appearance would suggest that sometimes the star simply 
spins about its symmetry axis without precessing. When 
the modulation was present, values in the range ~ 48° to 
~ 60° were found. Even if the spin-down were due entirely 
to a magnetic dipole, inequality (^^ shows that the phase 
residual can be extracted while neglecting magnetic torque 
variations. Then, if A"l> was small, equation (pfj) would give 
the wobble angle (up to a factor of tanx)- However, A<1? 
is not small, so the nonlinearised problem must be solved. 
This is not possible for a general value of x, but an example 
solution would be x = 59°, S = 80° for A* = 60° (equation 
p^). The key point is that 9 is not small, and the strains 
associated with such a large wobble angle would require a 
very high crustal breaking strain. Given that, to order of 
magnitude accuracy, the precession-induced strain is of or- 
der ^efluid, we have a strain of order 5 x 10~^. The breaking 
strain must be at least as large as this, exceeding even the 
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highest estimate of Ruderman (1992), weakening the free 
precession hypothesis. 

Now use equation ( |69| ) to extract the reference oblate- 
ness. Using a modulation period of 1009s we find eo ~ 
2.1 X 10~^A//_R. Using equation @ we find that the fiuid 
oblateness is 4.6x 10"^i?^/M. The inequality is violated 
(case I), even for high mass, small radius stars — the star is 
more nearly spherical than we would expect if its crust solid- 
ified at a rotation rate as high as 2.14ms. This discrepancy 
is resolved only if a significant portion of the star's total 
moment of inertia participates in the free precession. Using 
equation (^) we see that we need 7o = 0.337star for eo to be 
as large as efluid. 

Now we will test the gravitational wave spin-down hy- 
pothesis. Balancing the rate of loss of angular m omentum 
of the star against the gravitational fiux we find (Cutler & 
Jones 2|000|): 



32G^5(A7dr 

5c^ IstaT 



The quantity Aid can 
A7d//o = P/Pip, so that 



'(cos^6l + 16sin^6') (76) 
be extracted from observations: 



32G 5 

5c5 



P 

Wp 



-'o 



sin^ 6l(cos^ 6* + 16 sin^ 9) (77) 



Parameterising: 



/ 

20 fiHz/ day 



= 1.5 sin^ 0(cos^ ( 



16 sin 



/2.i4y / nr 

{~r) [pz 



lo 1 7sta 

0.33 



KhARl- (78) 



For Q \ and Jo = 0.33/star this is an order of magni- 
tude larger than the observed value. However, if we were 
to disregard the 6 ~ 1 derived from the phase modulation, 
we could use the above equation to calculate the wobble 
angle required for gravitational waves to provide angular 
momentum balance. We find Q — 0.44 radians = 25°. This 
corresponds to a crustal strain of order Stfluid ~ 2 x 10~^. 
The breaking strain must be at least as large as this. This 
is an uncor nfortably large v alue, exceeding even the highest 
estimates (Ruderman 1992). If the whole star participates 
in the free precession, so that Jo = Jstar, the wobble angle 
necessary to balance the spin-down is 9°, corresponding to a 
slightly more plausible breaking strain of 8 x 10""^. Clearly, 
if the free precession and gravitational wave spin-down hy- 
potheses are correct, our model points towards most or all 
of the star participating in the precession. 

One of the main pieces of evidence Middleditch et al. 
cite to support the gravitational wave spin-down hypothesis 
is that the frequency derivative / correlates rather well with 
the inverse square of the long modulation period Pfp. This 
is illustrated by figure 5 of Middleditch et al. (2000b), where 
/ varies by a factor of ~ 4, while Pfp varies by a factor of 2. 
As can be seen from equation (^) , this correlation is indeed 
predicted by this model, hut only for a fixed wobble angle. 
For the correlation to be perfect, the wobble angle would 
have to remain exactly constant. Referring to their figure, 
we see that with the exception of a single data point, all the 
error bars are consistent with the wobble angle remaining 
constant. The maximum fractional variation in the angular 
function sin^ d{cos^ d + lQs'm^ 9) consistent with / remaining 



within the error bars is approximately 10%, corresponding 
to a variation away from = 9° of less than 1°, for instance. 
It is not easy to see how the (unknown) process that leads 
to a factor of two variation in the size of the deformation 
AJd could preserve the wobble angle to better than 10%. 

To sum up, the large phase residual of around 60° can- 
not be accounted for: In such a fast spinning star it would 
correspond to a crustal breaking strain of around 5 x 10~^, 
an implausibly large value. Setting this difficulty aside, in 
order to satisfy inequality (^) we found that at least one 
third of the star had to participate in the free precession. In 
order for this motion to provide the necessary gravitational 
spin-down torque an even larger portion of the star must 
precess. If the whole star precesses, the gravitational wave 
spin-down hypothesis requires a crustal breaking strain at 
least as large as 8 x 10"'^. If a smaller portion precesses, an 
even larger breaking strain is required. 

7.3 Her X-1 

Her X-1 is a 1.24s X-ray pulsar in a 1.7d binary orbit about 
a main-sequence star. A third periodicity of 35d has been 
measured. Some authors have a ttributed this to forced pre- 
cession of the accretion disk ( Petterson 1975| ; Gerend & 
Boynton 1976; Bisnovatyi-Kogan et al. 1990), others to fre e 
precession of the star ( |Brecher 1972 ; Shakura et al. 1998 K 
while others argue for a com bination of both (Trumper et 
al. 198C; Ketsaris et al. 2000). If the 35d periodicity is due 
to free precession the wobb le angle must be large, of order 
unity ( ^hakura et al. 1998 ). 

Equation (|69|) gives a reference oblateness of eo = 
4.1 X 10~^, while equation (^ gives the star's actual shape as 
Efluid = 1-4 X 10~^. These values are consistent with inequal- 
ity ( p9| ) (case II), so that the observations are consistent 
with a crust-only precession with zero superfluid pinning. 
The pinned superfluid component can make up no more than 
about 10~* of the total moment of inertia. 



7.4 Summary and discussion of observations 

We will now comment upon our findings, which are sum- 
marised in table |l|. The key points are as follows: 

• Only one observation (Crab, Lyne et al. 1988) did not 
fit our free precession model for any choice of parameters — 
the calculated reference shape was too nearly spherical to 
have been formed by solidification earlier in the Crab's life 
when the star spun more rapidly than it does today. 

• Only one observation (Crab, Cadez & Galicic 1996a) 
required superfluid pinning for its explanation. Assuming 
crust-only precession, this pinned component must make up 
a fraction 10~® of the total stellar moment of inertia. 

• The remaining five observations could be made to fit our 
model with the pinned superfluid component set to zero, so 
that Coulomb forces provided the deformation in the mo- 
ment of inertia tensor. Two observations required at least 
part of the interior fluid to participate in the free preces- 
sion. Speciflcally, Vela required at least 6% of the total stel- 
lar moment of inertia to participate, while the SN 1987A 
remnant required at least 33%. The remaining observations 
were consistent with a crust-only free precession. The max- 
imum amount of pinned superfluid consistent with the ob- 
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Table 1. This table summarises stellar parameters calculated from the spin and (proposed) free precession periods. 
The quantity e^g is simply the ratio of these, P/Pfp. The actual oblatoness efluid is calculated using equation (^). 
The quantity €o is the reference oblateness as calculated when only the crust participates in the free precession, 
and there is no pinned superfluid (equation |69[ ). The quantity /sp/^star is the fraction of the stellar moment of 
inertia made up of pinned superfluid, assuming that only the crust participates in the free precession, and that the 
reference oblateness is zero (equation |ri| ). Key to references: (1) Lyne et al. (1988); (2) Cadez &; Galicic (1996a); 
(3) Deshpande & McCuUoch (1996); (4) Cordes (1993); (5) Stairs, Shemar & Lyne (2000); (6) Middleditch et al. 
(2000a); (7) Trumper et al. (1986). 



Object 


Reference 


Ecff 




Efluid 




£0 








B0531+21 (Crab) 


1 


6.4 X 10" 


-10 


2 X 10-* 




6 X lO"'^ 




9.6 X 10" 


12 


B0531+21 (Crab) 


2 


5.5 X 10- 


-4 


1.9 X 10" 


4 


5.5 X 10- 


1 


8.3 X 10- 


6 


B0833-45 (Vela) 


3 


6.2 X 10- 


-9 


2.7 X 10~ 


5 


6.2 X 10" 


6 


9.3 X 10- 


11 


B1642-03 


4 


4.5 X 10" 


-9 


1.3 X 10- 


6 


4.5 X 10" 


6 


6.8 X 10" 


11 


B1828-11 


5 


4.7 X 10" 


-9 


1.3 X 10" 


6 


4.7 X 10" 


6 


7.1 X 10- 


11 


SN 1987A remnant 


6 


2.1 X 10" 


-6 


4.6 X 10" 


2 


2.1 X 10" 


3 


3.2 X 10- 


8 


Her X-1 


7 


4.1 X 10" 


-7 


1.4 X 10" 


7 


4.1 X 10" 


4 


6.2 X 10" 


9 



servations is typically ~ 10^ of the stars' total moment of 
inertia. 

• We were able to explain two otherwise puzzling features 
of the Stairs et al. (2000) observation of PSR B1828-11 by 
assuming the dipole axis to lie nearly orthogonally to the 
deformation axis (x close to 7r/2). Specifically, we were able 
to account for the presence of a strong 504d periodicity in the 
1009d modulation, and the apparent discrepancy between 
the wobble angle as derived using the amplitude modulation 
and using the phase residuals. 

As stated, two observations require more than just the 
crust to follow the free precession. However, we wish to go 
further, and ask the following: Are the observations consis- 
tent with part (or all) of the core liquid in all the neutron 
stars following the free precession? Taking the extreme case 
of total crust-core coupling, we set 7o equal to the total stel- 
lar moment of inertia in equations (^8|) and ( |70| ) . This simply 
serves to increase the values of eo and /sf //star in the table 
by a factor of /star //crust, i.e. by a factor of ^ GlMi ^/R^ 
(equation 

In this case the qualitative conclusions that can be 
drawn are not very different from before. The Lyne et al. 
observation of the Crab is still inconsistent with the free 
precession model (although the reference oblateness is only 
a factor of 5 smaller than the fluid oblateness). The Cadez & 
Galicic observation of the Crab still requires a pinned super- 
fluid, now making up a fraction IQ-"^ of the total moment on 
inertia. As stated previously, the SN 1987A observation can 
only be explained assuming a large fluid component in the 
free precession. All the remaining observations continue to 
fit our model — the only change is that the maximum pinned 
superfiuid component allowed is increased to values typically 
of order IQ-*, still much smaller than the IQ-'^ predicted by 
theory. 

Another important question remains to be answered: 
Why are these stars precessing in the first place? What ex- 
cited this motion? The two main candidates for isolated stars 
are glitches and electromagnetic torques. A glitch could ex- 
cite free precession by occurring in a non-axisymmetric way, 
suddenly shifting the principal axis of the mome nt of inertia 



tens or w hile preserving the angular momentum ( Link et al _ 
199 ^). j Vlternatively, the spin-down electromagnetic torque 
can amplify the wobble angle of an already precessing star 



oldreich 1970). T here may ex ist also an 'anomalous' elec- 
tromagnetic torque (Jones 198S), which does not contribute 
to the spin-down, but will cause an initially non-precessing 
star to precess. Both glitches and electromagnetic torques 
are important for young stars, consistent with the rather 
young ages of the precession candidates. (All the isolated 
stars in table |l| have spin-down ages P/2p ^ 10^ years, 
with the exception of B1642-03, which has a spin-down age 
of 3 X lO'' years) . In the case of Her X-1, t he accretion torqu e 
is the obvious source of wobble excitation (Lamb et al. 1975). 

If electromagnetic torques are responsible, then we 
would expect almost all sufficiently young pulsars to display 
signs of free precession. In contrast, if glitching is respon- 
sible, we would expect to observe free precession in those 
stars which happened to have glitched recently, so that the 
precession has not yet been damped away. Only a fraction 
of the young pulsars show signs of free precession, favouring 
the glitch hypothesis. 



8 CONCLUSIONS 

As stated at the outset, our goal was to extract as much in- 
formation as possible from the handful of pulsar candidates 
for free precession. To this end we built a free precession 
model capturing (we hope) the main features of the prob- 
lem. Our model contained two coupling mechanisms between 
the crust and neutron fluid core. One was inertial coupling, 
where the fluid core effectively 'pushed' on the surrounding 
crust. The other was a frictional coupling due to the scat- 
tering of electrons off neutron vortices. We argued that even 
when these effects are taken into account, the neutron fluid 
core is not expected to follow the precession of the crust. 

When we compared our model against the observations 
we found that the wobble angle of the candidates was typi- 
cally small, less than a degree or so. Of greater interest was 
the information that might be extracted from the ratio of 
the spin and free precession periods, ecS- This ratio depends 
on the details of the star's structure, speciflcally on the ge- 
ological history of the crust (parameterised by eo), on the 
amount of pinned superfluid (/sf), and on the portion of the 
star that participates in the free precession (/o). It was not 
possible to examine the three parameters independently, as 
they all contribute to the observed ecS (see equation Bll) . 
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However, by assuming reasonable values for crustal 
strength (6) and breaking strain (ubroak), a few general con- 
clusions could be drawn. Firstly, if only the crust partici- 
pates in the free precession, superfluid pinning need be in- 
voked to explain only one of the observations. The others 
were consistent with no pinning at all, with an upper bound 
on IsF typically of 10"^" of the total stellar moment of iner- 
tia. In the (in our opinion less likely) case where the whole 
star precesses, the maximum pinned component was found 
to be typically 10~* of /star. Both these values are many 
orders of magnitude smaller than predicted by some glitch 
theories. Clearly, if the observations considered here really 
do represent free precession, superfluid pinning theory, at 
least as it affects free precession, is in radical need of re- 
working. 
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